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ABSTRACT 

We report the discovery of four groups of galaxies along the line of sight to the B1608-I-656 gravita- 
tional lens system. One group is at the redshift of the primary lensing galaxy (z = 0.631) and appears 
to have a low mass, with eight spectroscopically-confirmed members and an estimated velocity dis- 
persion of 150±60 km s~^. The three other groups are in the foreground of the lens. These groups 
contain ~10 confirmed members each, and are located at redshifts of 0.265, 0.426, and 0.520. Two of 
the three additional groups are centered roughly on the lens system, while the third is centered 
south of the lens. We investigate the effect of each of the four groups on the gravitational lensing 
potential of the B1608-f 656 system, with a particular focus on the implications for the value of Hq 
derived from this system. We find that each group provides an external convergence (next) of ~0.005- 
0.060, depending on the assumptions made in calculating the convergence. For lens systems where 
no additional observables that can break the mass-sheet degeneracy exist, the determination of Hq 
will be biased high by a factor of (1 — Hext)^^ if the external convergence is not properly included in 
the model. For the B1608-I-656 system, the stellar velocity dispersion of the lensing galaxy has been 
measured, thus breaking the mass-sheet degeneracy due to the group that is physically associated 
with the lens. The effect of the other groups along the line of sight can be folded into the overall 
uncertainties due to large-scale structure (LSS) along the line of sight. Because the B1608-I-656 system 
appears to lie along an overdense line of sight, the LSS will cause the measurement of Hq to be biased 
high for this system. The systematic bias introduced by LSS could be 5% or greater. Because LSS 
bias should be random along random lines of sight, averaging a large number of determinations of 
Hq from different lens systems should substantially reduce the uncertainties due to LSS on the global 
lens- based measurement of Hq. 

Subject headings: distance scale — galaxies: individual (CLASS B1608-f 656) — gravitational lensing 
- galaxies: groups 



1. INTRODUCTION 

One of the many possible uses of gravitational lens sys- 
tems is the measurement of Hq {ilefsdal 1964). However, 
this method has been limited by a lack of knowledge 
about the mass distribution in the lens, which leads to 
degeneracies between model parameters and Hq. Per- 
haps the most important degeneracy is that between the 
radial slope of the lensing mass profile and Hq. How- 
ever, strong limits can be placed on the mass slope in 
systems in which measureme nts of stellar velocity dis - 
persions can be made (e.g., iTreu fc Kooomaii^ I2004D . 
more than o ne component of the source is multiply- 
imag ed ('e.g-ICohn et al.ll2nf)H) . an Einstein ring is seen 
(e.g., lKocharie^'et'ldTl2001f) . or VLBI structure in the 

* BASED IN PART ON OBSERVATIONS MADE WITH THE 
NASA/ESA HUBBLE SPACE TELESCOPE, OBTAINED AT 
THE SPACE TELESCOPE SCIENCE INSTITUTE, WHICH IS 
OPERATED BY THE ASSOCIATION OF UNIVERSITIES FOR 
RESEARCH IN ASTRONOMY, INC., UNDER NASA CON- 
TRACT NAS 5-26555. THESE OBSERVATIONS ARE ASSO- 
CIATED WITH PROGRAM #00-10158. 
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lensed images c an be cl early mapped from one image 
to another fe.g.. lRusin et al. 2002 ). Another major de- 
generacy is that caused by an extended mass distribu- 
tion that is associated with the lens . This is the fa- 
mous "mass-sheet degeneracy" (e.g., iFalco et alllTflSl 
and can be caused by a cluster or group along the line 
of sight to the main lensing galaxy. The problem is 
that the standard lens observables, namely the locations 
and fluxes of the lensed images and the location of the 
lensing galaxy, do not indicate how much of the lens- 
ing mass surface density is due to the mass sheet. In 
fact, with the exception of very large sepa ration lenses 
(e.g. SPSS J1004+4112 and Q0957-H56 1: lOguri et alJ 
l200llWalsh. Carswell. fc WevmannI 11971 . it is difficult 
to know from the standard lensing observables whether 
or not there even exists an associated group or cluster. 
Thus, it is necessary to search for such structures by 
other means. We are conducting a survey of lenses in 
which time delays have been measured, with the aim of 
detecting groups or clusters which can affect the deter- 
mination of Hq from those systems. 
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The CLASS B1608+656 system lIMvers et al.lllQO.'iD is 
an excellent target for our survey. The redshifts of the 
lens and background source have be en measured to be 
z = 0.630 and z = 1.39, respectively l|Mvers et alJll995t 
iFassnacht et al .' 19961 It remains the only four-image 
lens system for which robust and high-precision measure- 
ments of all three independ ent time delays have been 
made (jFassnacht et al.ll2002j) . The lens has been sub- 
jected to intensive modeling, which incorporated infor- 
mation from the measured stellar velocity distribution 
and the Einstein ring shape. With an advanced modeling 
code, strong limits were placed on the slope of the mass 
density profile, yielding a measurement of Hq = 75lg km 
s~^ Mpc~^ (|Koopmans et all 12003). Furthermore, the 
mass models for this system requi re a relat ively large 
external shear of nearly 0.1 ^Koopmans et al.ti2003.'l . in- 
dicating the presence of nearby mass. In this paper, we 
provide evidence for a group of galaxies associated with 
the primary lensing galaxy and discuss the result that it 
and other structures along the line of sight have on the 
determination of from this system. 

Throughout this paper we assume r2j\/ — 0.3, JIa = 

0. 7, and, unless otherwise stated, we will express the 
Hubble Constant as Hq = lOOh kms"^ Mpc"\ 

2. OBSERVATIONS AND DATA REDUCTION 

We have conducted a spectroscopic survey of the 
B1608-I-656 field as part of our ongoing program to find 
compact groups of galaxies associated with gravitational 
lenses. The field wa s imaged in three bands, Gunn g, r, 
and i ijThuan & Cunn 1976). with the Palomar 60-Inch 
Telescope. Spectroscopic targets were selected based on 
their colors and distances from the lens system. The 
highest priority targets were those close to the lens sys- 
tem and with (r — i) colors in the range 0.45 to 0.65, 

1. e., close to those expected for early- type galaxies at the 
redshift of the lensing galaxy. All of the spectroscopic 
targets had r < 23. In addition to the high priority 
targets, several other galaxies were observed in order to 
pack efficiently the slitmasks that were used for the bulk 
of the spectroscopy. The spectroscopic observations were 
made with the Low R esolution Imaging Spectrograph 
(LRIS: lOke et al.llT995|l . in both longslit and multisht 
modes, and the Ech ellete Spectrograph and Imager (ESI; 
iSheinis et al."2002l on the W. M. Keck Telescopes. The 
spectroscopic data were reduced using scripts based on 
standard iraf^ tasks. Additional IDL scripts were used 
to process the ESI data. The spectroscopic and photo- 
metric data on all of the surveyed galaxies, along with 
the full details of the data reduction procedures, will be 
presented in a future paper dealing with the properties 
of the groups and the galaxies within them. 

3. MASS ALONG THE LINE OF SIGHT TO B1608+656 

The spectroscopic observations produced redshifts for 
97 galaxies in the field. The typical uncertainties in 
the redshifts, based on the scatter of redshifts calculated 
from individual lines in each spectrum, were Cz ~ 0.0003. 
The distribution of redshifts, in which several spikes are 

^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy 
under cooperative agreement with the National Science Founda- 
tion. 



seen, is shown in Figure^ We select group or cluster can- 
didates by looking for structures that are concentrated 
both spatially and in redshift sp ace, using the iterative 
method described bv lWilman et al. (2005). We investi- 
gated structures with at least five members in one red- 
shift bin, where the bin size was Az = 0.005. 

To characterize each group, we estimate its line-of-sight 
velocity dispersion ((T„) utilizing the ROSTAT package 
(Beers, Flynn, & Gebhardt 1990). This package avoids 
assumptions of gaussianity in the underlying velocity dis- 
tributions and includes bootstrap and jackknife error es- 
timation for the group redshifts and dispersions. Many of 
these estimators are significantly more resistant to non- 
gaussianity in the sample than the traditional Gaussian; 
large differences between the different scale measures can 
be a sign of deviations from gaussianity in the sample, 
either in shape or due to outliers. In particular, we 
use three estimators. The first, agap, is based on the 
gapper al gorithm, which is recommended for very small 
datasets (Beers et al."T990'). The error on this estimator 
is the jackknifed bi weight estimate. The second is the 
biweight scale estimate, abiwt , which is recommended for 
datasets with ~ 10 members. Its confidence interval is 
taken from the jackknifed gapper. Finally, we use cr gauss-, 
the Gaussian estimator, and its 1-sigma Student's t error. 
Consistency among the various dispersion estimates for 
a small redshift sample lends credence to the measure- 
ment. However, we note that velocity dispersions based 
on a small number of redshifts ca n be highly uncertain 
fe.g..lZa,b]udoff fc Mulcha,evllT99l. 

3.1. The Group Associated with B1608+656 

The most prominent spike in the redshift distribution 
consists of eight galaxies with redshifts of z 0.63. This 
spike includes the lensing galaxy at zi = 0.630. How- 
ever, we have not included the second lensing galaxy 
within the ring of imag es (G2; Kooum ans et alJ 120031: 
iSurpj fc Blandfordir2003ji . All indications are that G2 is 
merging with Gl, and therefore that this group consists 
of at least nine members. As one test of the likelihood 
that this redshift spike represents a real group, we have 
plotted the spatial distribution of the galaxies in the spike 
(Figure 2). We find that the galaxies are spatially con- 
centrated and centered roughly on the position of the 
lens. Seven of the eight galaxies in the redshift spike, 
or eight out of nine if G2 is included, are within a circle 
centered on the lens system of radius 2.'1, corresponding 
to 1 comoving Mpc at the redshift of the lens. Fig- 
ure 3 shows the redshift distribution in the region of the 
spike in terms of velocity offsets from the mean redshift 
of the group. The galaxies also have a tight distribution 
in velocity space, with all eight within ±300 km s^^ of 
the mean redshift. Thus, we conclude that these galaxies 
represent a group, hereafter called Group 1, associated 
with the lensing galaxy. 

For Group 1, we find that the various cr„ estimators 
obtained from the ROSTAT package are all consistent. 
This result strongly suggests that, despite the small num- 
ber of redshifts in the group, the velocity distribution is 
well described by a Gaussian, and the resulting veloc- 
ity dispersion is robust. In Table we list the median 
redshift and the three estimates of the line-of-sight veloc- 
ity dispersion {<Jgap, ouwt, and u gauss) with their associ- 
ated errors (in km s~^). Based on the number of group 
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Fig. 1. — Histogram showing the distribution of the 97 non- 
stellar redshifts obtained in the field of B1608-I-656. The width 
of the bins is Az = 0.005. The most prominent spike, with eight 
galaxies, is at the redshift of the lensing galaxy (z = 0.63). 



Fig. 3. — Histogram showing the velocity distribution of the 
galaxies in and surrounding the redshift spike at z = 0.63. The bin 
widths are 200 km s~^, approximately twice the uncertainties in 
determining the redshifts. The curve is a Gaussian with cr =150 km 
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Fig. 2. — Spatial distribution of the galaxies in Group 1. The 
field of view is lO'x 11', with the axes labeled in terms of offsets 
from the B1608-I-656 lens system in arcseconds. The dots represent 
the positions of galaxies with r < 23, while the open circles mark 
the galaxies for which redshifts have been obtained. The open 
boxes mark the galaxies in the group. The large dashed circle has 
a radius of 1 comoving Mpc at the redshift of the lensing 

galaxy. 



members, we use cr„ = a gap for Group 1, since a-gap is 
the more appropri ate estimator for very small data sets 
ijBeers et al.lll99(H) . Thus, we obtain = 150 ± 60 km 



3.2. Other Mass Concentrations Along the Line of Sight 

There are other spikes in the distribution of measured 
redshifts shown in Figure ^ In addition to the group at 
the redshift of the lens, there are three other group candi- 
dates in the observed distribution satisfying the redshift 
and spatial concentration criteria, with mean redshifts of 



< z 0.265 (Group 2), < z 0.426 (Group 3), and 

< z >= 0.520 (Group 4). Each group has a substantial 
number of confirmed members, with sizes of nine, seven, 
and 14 galaxies, respectively. The properties of these 
groups are given in Tabled For each group, we find that 
the various estimators are all consistent. Based on the 
number of confirmed members in each group, we estimate 
the group velocity dispersions with agap for Groups 2 and 
3, and abi^t for Group 4. Figures ^ and El show the spa- 
tial and velocity distributions, respectively, for the three 
additional groups detected in this field. We explore the 
effects that these groups may have on the determination 
of Hq from B1608-I-656 in the following sections. 



4. DISCUSSION 
4.1. A Cluster at z = 0.52? 

Of the three additional candidate groups, the one at 
z = 0.520 stands out. Its velocity dispersion of ~900 km 
s^^ implies that this structure is a cluster of galaxies 
rather than a group. If it is in fact a real cluster, it 
will have a significant effect on the measured value of 
Hq obtained from the B 1608-1-656 lens system (see next 
section). Therefore, we consider evidence that pertains 
to the reality of the cluster. 

There are two main arguments in favor of the redshift 
spike's representing a real cluster. The first is the con- 
centration of the structure in both velocity space and in 
projection on the sky. Nine of the 14 spectroscopically- 
confirmed galaxies in the structure are within ± 1000km 
s~^ of the mean redshift, while 13 of 14 lie within a pro- 
jected distance of 1 comoving Mpc of the median 
position. The second is that nine of the 14 galaxies in 
the group lie along a tight sequence on a color-magnitude 
diagram, at (r — ?) ~ 0.45 (Figure EI)- 
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Fig. 4. — Spatial distributions of the galaxies in the additional three group candidates in the B1608+656 field, represented as in Figurel^ 
The open circles represent galaxies for which redshifts have been obtained, while the open boxes represent confirmed group members. The 
dashed circle in each plot has a radius of lh~^ comoving Mpc at the redshift of the group and is centered at the median position of the 
confirmed group members, (a) Group 2 at 2 ~ 0.27. (b) Group 3 at 2: 0.43. (c) Group 4 at 2 ~ 0.52. 



TABLE 1 
Group Properties 



^ rr rr ^ PA ^ fi ^ PA t) 

^gap ^biwt ^ gauss ^med ^ ^med ^Iw ^ ^Iw 

Group < 2 > Di^/Ds Ngais (km s"-*-) (km s"-') (km s"-*-) (arcsec) (°) (arcsec) (°) 



1 


0.6313 


0.45 




150±60 


130 ± 100 


150t^o 


25 


30 


8 


-5 


2 


0.2651 


0.74 


9 


320±100 


230 ± 270 


320ilgo 


69 


29 


18 


6 


3 


0.4263 


0.66 


7 


270±110 


260 ± 120 


270 ± 110 


63 


-166 


46 


-155 


4 


0.5202 


0.53 


14 


920±120 


930 ± 140 


890l?« 


13 


26 


6 


-47 


4ad 


0.5163 


0.53 


7 


410±160 


100±320 


430±110 


24 


47 


46 


6 


4bd 


0.5241 


0.53 


7 


360±120 


350±120 


350±120 


44 


-150 


52 


-161 



^Median group position expressed as an offset from the lens. The PA is measured north through cast. 

^The luminosity-weighted group position expressed as an offset from the lens. The PA is measured north through east. 

*^This number does not include the second galaxy (G2) inside the Einstein ring of the lens system because there is no 
measured redshift for G2. However, G2 appears to be merging with the primary lens galaxy and thus appears to be a 
member of the group. 

^Groups 4a and 4b arc subsets of Group 4. See §4.1. 



The arguments against the reahty of the cluster are as 
foUows. First, an examination of deep Advanced Cam- 
era for Surveys (ACS) images of this field (Program GO- 
10158; PI: Fassnacht) shows no obvious overdensity such 
as might be expected from a cluster. Furthermore, of 
the 11 galaxies that are members of Group 4 and also 
lie within the field of view of the ACS imaging, none ap- 
pears to be a central bright early- type galaxy. In fact, 
nine out of the 11 group galaxies in the ACS imaging 
appear to be spirals. This is surprising given the red col- 
ors of these galaxies, but we will leave the discussion of 
this point to the foUowup paper in which we will focus 
on the properties of the groups and the galaxies within 
them. The second point arguing against the cluster hy- 
pothesis is that the group velocities (Figure 5c) are not 
normally distributed as one would expect from a relaxed 
cluster. With the large velocity dispersion implied by 
the measured redshifts, it is perhaps not unlikely that a 
random set of 14 galaxies does not produce a nice Gaus- 
sian shape. Further spectroscopy of the field, especially 
pushing deeper than the r ~ 23 limit used in the previous 
observations, may fill in the velocity structure and pro- 
duce a velocity distribution that is closer to Gaussian. 
On the other hand, the last slitmask in our observing 
program had the targets re-prioritized to preferentially 



find galaxies in the cluster. Only one additional member 
was found, once again indicating that there is not the 
large overdensity of galaxies typically found in a cluster. 
Third, there was no extended X-ray emission detected in 
this field in a recent 30 k sec Chandra X-ray Obs ervatory 
observation of this field l|Dai fc Kochanel3l200!Tt) . If the 
cluster is as massive as suggested by its velocity disper- 
sion and is virialized, it should have been detected by 
the Chandra observations, which had a 3-tT upper limit 
for detection t hat correspond to a velo city dispersion of 
-500 km s"^ l|Dai fc Kochanekl |2005|) . However, it is 
possible that the cluster is recently formed or is, in fact, 
a pair of merging groups. In these situations, it may be 
severely under luminous at X-ray wavelengths. Finally, 
a cluster this massive and centered as close to the lens 
as it appears to be would introduce a strong effect on 
the lensing signature. The external shear introduced by 
such a cluster would b e much la rger than required by 
the lens model (Koopmans et al. l i200.3.) . Furthermore, it 
is possible that the surface mass density would be high 
enough at the position of the lens to form yet another 
image of the lensed source, i.e., the cluster would have 
to be treated in the strong lensing regime. 

We feel that the arguments against the reality of the 
cluster are stronger than those in favor. The structure 
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Fig. 5. — Histograms showing the velocity distributions of the galaxies in the three other spatially-concentrated rcdshift spikes. The bin 
widths are 200 kms~^, approximately twice the uncertainties in determining the redshifts. 
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Fig. 6. — Color-magnitude diagram for the B1608+656 field. 
The dashed lines represent the approximate completeness limit of 
the imaging obtained at the Palomar 60-Inch Telescope. The boxed 
points belong to the z = 0.52 cluster candidate. The diagonal solid 
line represents the limit of the spectroscopy at r ~ 23. 



may instead be a pair of merging groups or some kind 
of filamentary structure. To explore the two-group hy- 
pothesis, we arbitrarily split the group into two parts 
by redshift. Group 4a is defined as the seven galaxies 
with redshifts smaller than the median Group 4 redshift, 
while Group 4b is comprised of the seven galaxies with 
redshifts larger than the median. Figure [7| shows the 
spatial and velocity distributions of the two groups. The 
velocity dispersions of the two groups are ~ 400 km s~^ 
and ~ 350 km s~^. If the two-group explanation is cor- 
rect, and these velocity dispersions are close to the true 
values, t hen it is not surp rising that the X-ray observa- 
tions of iDai fc Kochaiiela l|2005.1 did not detect diffuse 
X-ray emission. The spatial distribution of Group 4b is 
centered slightly to the southwest of the Group 4a cen- 
troid. However, the two centroids are completely con- 
sistent, given the uncertainties. Another possibility is 
that the Group 4 galaxies lie along a filament. We note 
that much of the following discussion treats each group 



as a collection of individual halos, each associated with a 
group galaxy. Therefore, whether the galaxy under dis- 
cussion is in a cluster, a group, or a filament is irrelevant. 

4.2. Effect on Gravitational Lensing 

We now consider the effect the groups associated with 
B1608-I-656 may have on the overall lensing gravitational 
potential, with a particular emphasis on the effect on Hq. 
While the angular separation of the lensed images in a 
gravitational lens system provides an accurate measure- 
ment of the projected lensing mass, it does not require 
that all of the mass be associated with the primary lens- 
ing galaxy. In fact, some of the mass can be contributed 
by other structures along the line of sight, such as the 
groups discussed above. The contribution of external 
mass is quantified through the convergence, Hext, that it 
causes at the location of the lensing galaxy. The conver- 
gence is just the scaled mass surface density: 



At: G Df D 



is 

where Sc represents the critical mass surface density re- 
quired for multiple images to form. As usual, ZJ^s, 
and Ds are the angular diameter distances between ob- 
server and lens, lens and source, and observer and source, 
respectively. The external convergence will lead to a 
value of i/o that is too high if the full lensing mass is 
improperly assigned solely to the primary lensing galaxy, 
i.e., 

where H^ -meas is the value obtained without properly in- 
cluding the external convergence in the lens model. We 
will use two methods to estimate the external conver- 
gence contributed by each of the groups along the line 
of sight to the B1608-I-656 system. Although neither of 
these methods may be entirely correct, the range of re- 
sults that they produce should be representative of the 
true group convergences. 

The first, and more traditional, method is to assume 
that the group can be approximated as a smooth mass 
distribution. The distribution, for simplicity, is usu- 
ally taken to be that produced by a singular isothermal 
sphere (SIS). In this case, the convergence contributed 
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by the group, calculated at the location of the lens, is 
Des 2-k(jI bsis 



K,SIS(fcent) 



29, 



where Ocent is the angular offset between the cen- 
ter of the group and the lens system. The "lens 
strength" of an isothermal distribution is defined as 6 = 
Ancr^Dgs/iDsC^), and for a singular isothermal sphere 
gives the Einstein ring radius in angular units. If Ocent is 
measured in arcminutes and is measured in km s~^, 
then 



Ocent \ 

arcmin / 



The convergences due to the cluster and group candi- 
dates, calculated using the SIS assumption, are given 
in Table |21 In each case the velocity dispersion used 
was that obtained from the gapper method. With the 
strong dependence of Kext on Ocent, it becomes imper- 
ative to locate the group centroid accurately, which is 
extremely difficult with fewer than several tens of con- 
firmed group members. For example, the centroid esti- 
mates listed in Tabled differ substantially depending on 
whether they were obtained by taking the median posi- 
tion or the luminosity-weighted position. Therefore, we 
also apply a second method for estimating the group con- 
vergence. 



The alternate method for computing the group conver- 
gence is to treat the group as a collection of individual 
galaxy halos, with no overall group halo. In other words, 
the mass sheet at the position of the lensing galaxy is 
composed of the overlapping halos of the other galaxies 
in the group. The use of this method is motivated by two 
considerations. First, by ignoring any shared group halo, 
this method explores what is probably an extreme case. 



especially when assuming that the galaxy halos may be 
truncated (see below) . This extreme case should provide 
a fairly robust lower limit to the estimate of the group 
convergence. Second, these moderate-redshift groups 
may still be in the early phases of formation and thus 
the galaxies may not yet have lost a significant fraction 
of their individual h alos to the shared group ha lo. Simu- 
lations performed bv lKeeton fc Zabludofj (2004) indicate 
that similar results are obtained whether the groups are 
treated as a single shared halo or a collection of individual 
halos. Furthermore, the input cosmological parameters 
(e.g., Hq) are recovered accurately when the group con- 
tribution is treated as a collection of galaxy halos. This 
accuracy is obtained even though they make the simplify- 
ing assumption that the group galaxies are circular when, 
in fact, their simulated galaxy mass distributions were el- 
liptica l. We follow the approach of iKeeton &: ZabludofJ 
^|20n4^ and assign group galaxy masses (expressed in 
terms of their lens strengths, bi) based on their optical 
magnitudes (rrii) via the relationship 

where b fid and m fid are the lens strength and magnitude 
for a fiducial galaxy within the group. In each case, we 
take the fiducial galaxy to be the brightest confirmed 
group member. If each of the grou p members i s treated 
as a si ngular isothermal sphere, as in lKeeton fc Z abludo3 
l)2004|) . then the total convergence at the position of the 
lens is just the sum of the individual convergences at that 
position: 

E_ \ " bi 

i z * 

where Oi is the distance between galaxy i and the lens. 
The results for each of the four groups are given in Ta- 
ble |21 while notes on the individual groups are given be- 
low. 

Another effect to consider is that due to the truncation 
of the dark matter halos of the group galaxies. In the cal- 
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TABLE 2 

Convergences due to Groups 



Group 








f^trunc 


AT d 


1 


0.0056 


0.018 


0.025 


0.012 


1 


2 


0.016 


0.061 


0.013-0.026 


0.0040-0.0082 


1 


3 


0.011 


0.015 


0.014-0.028 


0.0064-0.013 


2 


4 






0.026-0.053 


0.015-0.031 


3 


4a 


0.054 


0.028 








4b 


0.023 


0.019 









^Convergence calculated with the group represented as a SIS and the group 
centroid represented by the median galaxy position. 

^Convergence calculated with the group represented as a SIS and the group 
centroid represented by the luminosity-weighted mean. 

^Rangc corresponds to a range of stellar velocity dispersions for the fiducial 
(brightest) galaxy in Groups 2, 3, and 4. The velocity dispersions range from 
140-200 km s~^ 

'^Number of galaxies contributing to /^trunc- 



culation of Kind, all of the galaxy halos are assumed to be 
larger than the separation between the galaxies and the 
primary lens system. However, w eak lensing studies by 
iHoekstra. Yee. fc Gladder^ l)2004f) have suggested that 
galaxy halos have a truncation radius of ~ 200h~^ kpc. 
If this truncation radius is real and is typical, then the 
convergence due to galaxies located farther from the lens 
than this will drop off faster than the (1/9) assumed in 
the calculation of Kind- To explore the possible size of 
this effect, we recalculate the external convergence con- 
tributed by each group, summing only the contributions 
from the galaxies that lie within a projected distance 
of 2QQh~^ comoving kpc from the lens. The results are 
given as Ktmnc in Table [3 This is an extreme approach, 
but it should give an approximate lower limit to the con- 
vergence caused by each group. 

4.2.1. Convergence Due to Group 1 

For this group we use F814W magnitudes measured 
from Hubble Space Telescope observations of the field. 
These observations consist of the deep ACS images men- 
tioned above, as well as Wide-Field Planetary Cam- 
era 2 images of the same field (Program GO-6555: PI 
Schechter). The fiducial galaxy is the primary lensing 
galaxy (Gl) which has a F814W magnitude of 18.2 and 
a lens strength of 616O8 = 0'.'83 (Kooomans ct al. 2003). 
The group galaxy most distant from the lens system is 
not covered by the HST images so we do not include it 
in the calculation of Kind- However, it is far enough away 
(2.'7) that its contribution to the overall convergence is 
negligible. 

4.2.2. Convergence Due to Croup 4 

As discussed in §4.1, we do not beheve that Group 4 
is a real cluster. We note that an assumption that the 
cluster is real yields extremely large convergences at the 
position of the B1608+656 system. The SIS approxima- 
tion leads to ksis 6.5/0. For an isothermal sphere, an- 
other image of the background object will be produced 
when K > 0.5. This condition is satisfied for the fore- 
ground cluster whether the cluster centroid is estimated 
using the luminosity-weighted position or the straight 
median position. Using these centroids leads to k = 1.1 
or K — 0.51, respectively, at the position of the lens. We 



have not seen evidence for an obvious lensed counterpart 
in deep 5 GHz or 8.5 GHz radio maps, although there 
are other compact so urces in the field (e.g., object 2; 
iFassnacht et al.lll999|) . Furthermore, such large conver- 
gences, whether from a SIS or some other mass distribu- 
tion would almost certainly lead to an image separation 
in the main lens larger than the 2'.'1 that is observed. 
Once again, these arguments indicate that Group 4 is 
not a real cluster. Because the evidence against the re- 
ality of the cluster is strong, we instead apply the SIS 
approximation to our arbitrarily selected subgroups, 4a 
and 4b. For these groups, the convergences (Tabled) are 
much more reasonable 

The calculation of Kind for Group 4 does not depend 
on whether the group is a single cluster or a pair of merg- 
ing groups or a filament. For this group, the galaxy 
lens strengths are estimated from the r-band magnitudes 
from the ground-based imaging. The fiducial galaxy, 
which is the brightest confirmed group member, has 
r = 20.0. We will discuss different estimates of the lens 
strength of the fiducial galaxy below. We note here, how- 
ever, that a galaxy at this redshift will have 6 '--^ 1 if its 
stellar velocity dispersion is '^260 km s~^. 

4.2.3. Convergence Due to Croups 2 and 3 

The convergences calculated for Groups 2 and 3 using 
the SIS method are larger than that of Group 1 due to 
their larger velocity dispersions. For Group 3, the two 
SIS estimates of the convergences are similar because the 
group is compact and located approximately 1' from the 
lensing galaxy. In contrast, because Group 2 is roughly 
centered on the lens system, small displacements in the 
centroid can lead to large changes in ksis- For both 
Group 2 and Group 3, the fiducial galaxies used in the 
Kind and Ktrunc mcthods are once again the brightest 
galaxies in their respective groups. In each case the fidu- 
cial galaxy has r = 20.0. Galaxies at the redshifts of 
Groups 2 and 3 have 5 ~ 1 for stellar velocity disper- 
sions of ^220 km s~^ and '^230 km s^^, respectively. 

4.2.4. Total Convergences of Croups 2, 3, and 4 

In order to compute the convergences for Groups 2- 
4, we need to assign a value for bfid for each group. 
We do this by assuming that the fiducial galaxy in each 
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group has a velocity dispersion of 200 km s^^. This is 
sUghtly smaller than the expected dark matter velocity 
dispersion for a L* galaxy of ~ 225 km s~^ (e.g., 

^ochanck 1996). Because the fiducial galaxy in each case 
is the brightest one in its respective group, this assump- 
tion should be reasonable. The resulting lens strengths 
are 0.85, 0.76, and 0.61 for Groups 2, 3, and 4, respec- 
tively. The resulting convergences are given as the upper 
end of the ranges given for Kind and Ktrunc in Table El 
These are within a factor of a few of the convergences es- 
timated from assuming that the groups were isothermal 
spheres. 

One factor that could lessen the total convergence 
arises from the conversion of the group galaxy luminosi- 
ties to masses. For example, in the ACS images, several 
of the confirmed members of the groups appear to be 
late-type galaxies. In this case the estimated masses are 
probably too high since spir als have lower mass-to-light 
ratios than ellipticals (e.g., IBahcall. Lubin. fc DormanI 
Il995jl . To explore this effect, we assume that the fiducial 
galaxy for Groups 2, 3, and 4 has a velocity dispersion 
of 140 km s""'^. The resulting lens strengths are ~50% as 
large as those obtained above. The corresponding con- 
vergences are given as the lower ends of the ranges listed 
in the Kind and Ktrunc columns of Table |2| 

In contrast, it is almost certain that not all of the group 
members have been identified. An incomplete census of 
the group will lead to an underestimate of the total con- 
vergence. We note that the deep ACS imaging of this 
field reveals at least 10 galaxies within 10" of the lens, 
corresponding to ~ 40 — 90 h^^ comoving kpc at red- 
shifts between 0.3 and 1.0. Only one of these galaxies 
has a measured redshift {z = 0.6087). Overall, a reason- 
able estimate is that true convergence from each group 
falls between the lower end of the Ktmnc range and the 
largest of the other estimates of Kext. 

4.3. Stellar Velocity Dispersion of the Lensing Galaxy 

Until now we have ignored another datum which can 
be used to break the mass-sheet degeneracy, namely the 
measured stellar velocity dispersion of the lensing galaxy. 
This measurement provides an estimate of the enclosed 
mass at a smaller radius than that probed by the lens- 
ing signature. Thus, the mass measurement from lens- 
ing can be combined with that from stellar dynamics 
to provide an effe ctive mass density slope in the lensing 
galax y (e.g.. lKooDmans fc T rcu 2002: iTreu fc Kooomana 
I2002D . The result of adding an external mass sheet that is 
physically associated with the lensing galaxy is to flatten 
the overall mass density profile, which leads to a lower 
value of Hq for given time delays. A measurement of the 
stellar velocity dispersion will reflect this flattening, i.e., 
the velocity dispersion will not be as high as would have 
been predicted from the lensing mass and an assumption 
of a steeper mass density profile. For a density profile 
that is close to isothermal (i.e., p oc r^^), essentially the 
same value of Hq should be derived whether the system 
is modeled as a single galaxy with the effective mass den- 
sity slope measured from lensing plus dynamics or as a 
galaxy with a steeper density profile but also including 
an external mass sheet (Koopmans 2004). For a high- 
accuracy measurement of Hq , it is thus crucial to obtain 
high precision measurements of either the external con- 
vergence or the stellar velocity dispersion of the lensing 



galaxy. In the case of B1608-f 656, the measured velocity 
dispersion is 247 ± 35 km s^^, indicating that the mass 
distribution in the lens (including any contribution from 
the group) is close to isothermal ( Kooomans et al. 2003) . 
The uncertainties in the velocity dispersion are, in fact, 
the largest source of error in the current model. 

4.4. Effect of Large-scale Structure 

Finally, we consider the effect of the mass in large-scale 
structures (LSS) along the line of sight to the background 
source. In an ideal situation, it would be possible to ac- 
count for all mass along the line of sight and to trace 
the rays from the background source through the distri- 
bution. In practice, even with deep space-based imag- 
ing, this is impractical. The uncertainties arising from 
photometric redshifts, the conversion of light to mass, 
etc., would far exceed the expected size of the effect. 
Therefore, it is necessary to examine the effect of LSS 
in a statistical manner. Analytic calculations have in- 
dicated that large-scale structure should affect the time 
delay s for a lens system, and hence Hq, by a few per- 
cent (jSeliaklll994f) or perhaps as much as 10 %, depend- 
ing o n the redshift of the background source l)Bar-Kanal 
IT99I . The effect of LSS can be to either increase or 
decrease the value of Hq because voids along the line 
of sight effectively act as areas of negative density when 
compared to the mean density of the Universe. Although 
this effect is a systematic one for any given lens sys tem, it 
shoul d be random for random lines of sight (e.g.. ISeliaS 
I1994|) . Therefore, it should be possible to significantly 
reduce the uncertainty due to large-scale structures on 
the global measurement of Hq by averaging the values 
obtained from many lens systems. We note that lenses 
may lie along lines of sight that are biased (e.g., have 
more line-of-sight structure) and therefore that the effect 
of large-scale structure can not be eliminated completely 
by averaging lens-based measurements of Hq. We will 
investigate this question in a future paper (Fassnacht et 
al., in prep). 

4.4.1. The Effect on Hq 

Finally, we arrive at the total effect on the value of Hq 
derived from the B1608-I-656 system. The group physi- 
cally associated with the lens system. Group 1, provides 
an external convergence of ~ 0.01 — 0.03 (Table [^J. As 
discussed above, however, the degeneracy associated with 
Group 1 is broken by the measurement of the stellar ve- 
locity dispersion of the lensing galaxy. Therefore, the 
effect of Group 1 on the value of Hq derived from this sys- 
tem gets folded into the determination of the radial mass 
slope in the lensing galaxy, and is thus already included 
in the published uncertainties on Hq from B1608+656. 

The other groups along the line of sight also provide 
convergences of a few percent. These groups are almost 
certainly typical of the clumpy structure that can be 
found along any line of sight. None of them provide 
an extraordinarily large convergence, unless Group 4 is 
a real cluster. Therefore, one would be inclined to in- 
corporate the effects of these groups into the effects of 
LSS and conclude that these groups will contribute to 
the random shift of a few percent in the value of Hq de- 
termined from this system. However, as a result of our 
investigation of lens fields (Fassnacht et al. , in prep) we 
have determined that the B1608-f 656 system lies along 
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a line of sight that is overdense compared to typical lines 
of sight. Therefore, the expected effect of properly incor- 
porating the effects of LSS on this system is to reduce the 
value of Hq from the previously measured value. If the 
majority of the extra mass along this line of sight is being 
contributed by Groups 1-3, the size of the effect could be 
5% or more, given the convergences in Table |21 This ef- 
fect, thus, could be comparable in size to the statistical 
uncertainties quoted for the Hq measurement from this 
system, and would reduce the central value to '^70 km 
s~^ Mpc~^ or lower. 

5. SUMMARY AND FUTURE WORK 

Our spectroscopic observations of the field containing 
the lens system B1608+656 have provided evidence for 
four groups of galaxies along the line of sight to the lens 
system. These groups should contribute external conver- 
gence at the location of the lens and, therefore, will affect 
the determination of the Hubble Constant obtained with 
this system. However, quantifying the amount of exter- 
nal convergence is difficult. If each group contains all 
of its mass in a overall smooth halo, it becomes imper- 
ative to measure accurately the halo mass and centroid. 
These measurements can be highly uncertain when made 
based on ^ 10 redshifts. We therefore follow a second ap- 
proach and treat each group as a collection of individual 
galaxy halos. To establish a firm lower limit to the con- 
vergence due to each group, we examine the convergence 
contributed only by galaxies within a projected distance 
of 200 h"^ comoving kpc from the lens. However, our cal- 
culations do not include contributions from other galax- 
ies within the truncation radius of the lens because their 
redshifts are unknown. Each of the groups contributes a 
convergence of a few percent. The overall effect on Hq 
is less than suggested by the sum of the convergences 
because the stellar velocity dispersion of the main lens- 
ing galaxy has been measured. This measurement breaks 
the mass-sheet degeneracy, at least due to Group 1, by 
determining the effective radial mass density slope in the 
lensing galaxy. The effects of the other groups can be 
folded into the overall uncertainties due to large-scale 
structure along the line of sight, since none of the groups 
appears to be extremely massive. However, because the 
line of sight to B1608-I-656 appears to be significantly 
overdense compared to typical lines of sight, the effect 
of LSS will be to bias the current measurement of Hq 
high, i.e., the true value of Hq from this system should 
be lower than the published value. The size of the effect 
of the LSS could be 5% or larger. We note that none 
of the newly discovered groups was obvious from opti- 
cal images of this system. Therefore, it is important to 
closely investigate the fields of lens systems in order to 
examine possible sources of bias in lens-based measure- 
ments of Hq, especially if the stellar velocity dispersion 
of the lensing galaxy is unknown. 



We are actively working to reduce the uncertainties in 
the determination of Hq from this lens system. The deep 
ACS imaging is being used as an input to new lensing 
code that can properly incorporate Einstein ring struc- 
ture. This modeling should reduce the mass-slope uncer- 
tainties substantially in the region of the lensed images. 
We will also use the ACS imaging to search for weak 
lensing signatures of mass concentrations along the line 
of sight to the lens system. The weak lensing will pro- 
vide further constraints on the amount of external con- 
vergence for B1608-I-656. Another approach to reducing 
the uncertainties in the mass slope would be to obtain 
higher sensitivity spectroscopy of the lensing galaxy. The 
spectroscopy would provide a more accurate determina- 
tion of the stellar velocity dispersion. Further optical 
spectroscopy of the galaxies in the field and deeper X-ray 
imaging would also add information on the external mass 
distribution. Finally, by averaging measurements of Hq 
obtained from many different lens systems, we should be 
able to reduce the uncertainties due to large-scale struc- 
ture and thus obtain a precise global measurement of Hq 
from lensing. 
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